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^ ' Abstract 



In this paper, we study the properties of a phase shp in a superfiuid Fermi gas near a Feshbach 
resonance. The phase shp can be generated by the phase imprinting method. Below the superfiuid 
transition temperature, it appears as a dip in the density profile, and becomes more pronounced 
when the temperature is lowered. Therefore the phase slip can provide a direct evidence of the 
superfiuid state. The condensation energy of the superfiuid state can be extracted from the density 
profile of the phase slip, due to the unitary properties of the Fermi gas near the resonance. The 
width of the phase slip is proportional to the square root of the difference between the transition 
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\Q , temperature and the temperature. The signature of the phase slip in the density profile becomes 
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A lot of effort has been concentrated on creating the superfluid state in a Fermi gas [l| . 
In recent years, the quantum degenerate regime was reached [2] and Feshbach resonance was 



observed in Fermi gases 



s. 



j . Asymmetric expansions of cold Fermi gases near the resonance 
were observed |^, |5|, as a result of the strong interaction between atoms 5]. The molecular 
condensate was observed when the scattering length was tuned to the positive side of the 
.escance fl Q. Conde„.a.o„ of Fe^ion pa.s wa, observed whe. the Fe.n. gas was fa. 
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tuned through the resonance [8|, |9|| . The measurement of the breathing-mode frequency was 
consistent with the predication of superfluid hydrodynamics, and the damping of the mode 
was weak at low temperatures !10|. Across the Feshbach resonance, a sudden change in the 



collective excitation spectrum was observed 



0. 



More recently, observation of the pairing 



gap in the radio- frequency excitation spectra of the Fermi gas was reported [l^- However, 
phase coherence of the wave-function, one of the most important properties of the superfluid 
state, has not been observed in these experiments so far. In this paper, we propose that 
phase coherence of the superfluid state can be observed in the phase-slip experiment. 

Phase coherence is one of the key characteristics of macroscopic quantum phenomena 
such as Bose-Einstein condensation, superfluidity, and superconductivity. In these systems, 
the phase symmetry of the wavefunction is broken. When two systems are coupled to- 
gether, the two order parameters interfere with each other. The Josephson effect appears 
when two superconductors with different phases are linked by a junction. In a system with 
two Bose-Einstein condensates, the interference pattern was observed upon release of the 
condensates (l3l| . In quasi-lD systems, the order parameter in the middle region vanishes 
when two sides of the system have different phases. This phenomenon is called phase slip. 
In quasi-lD superconductors, phase slip is the main reason for flnite conductivity |1J| . In 

lej. They 



Bose-Einstein condensates, phase slips have been observed in experiments [l5|, 
are also often called solitons because they have the same dynamics as solitons. 

In the superfluid state of the Fermi gas, the phase slip can be generated by the phase 
imprinting method. A pulsed laser can illuminate the system with different intensities on 
the left and right sides, and the pulse duration can be tuned so that for an atom the laser 
produces a phase difference of 7r/2 or — 7r/2 between the two sides. For a pair of atoms, 
the phase difference is twice as much, vr or — vr, between the two sides of the system. As in 
superconductors, the order parameter in the superfluid state describes the pairing of atoms. 
When there is a discontinuity in the phase of a pair wave-function, the order parameter must 



vanish at the point of discontinuity to generate a phase shp. Phase imprinting methods have 
been successfully applied to create vortices [l7|, |l8| and solitons [15|, |l6| in Bose-Einstein 
condensates. However, in the normal state of a Fermi gas, there is no phase coherence, and 
the phase imprinting does not produce any special effect. In the following, we study the 
properties of a phase slip in a superfluid Fermi gas near a Feshbach resonance. It is assumed 
that the phase step in the phase imprinting is applied along the axial direction of the trap 
and the phase slip is generated at the center of the trap. 

In the superfluid state, the free energy of a Fermi gas is lower than that in the normal 
state. This energy difference is defined as the condensation energy. The total free energy L 
is the sum of the condensation energy Lc and the free energy of the normal state Ljv, 

L = Lc + Ln. (1) 

According to the Ginzburg-Landau theory, the condensation energy can be expanded in 
terms of the order parameter which is proportional to the energy gap of the fermions. 
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(2) 



where a and (3 are coefficients. Below the transition temperature Tc, the sign of the coeffi- 
cient a is negative and the superfluid state has lower free energy than the normal state. The 
free energy of the normal state is generally a function of temperature, particle density, trap- 
ping energy, and interaction energy which is usually proportional to the scattering length. 
The total free energy is a minimum in the space of the particle density and the superfluid 
order parameter, 

^ = »' (^' 

.^^0. (4) 



where n(r) is the fermion particle density at location r. The derivative in equation Q is 
taken under the constraint of constant total-particle number. The equilibrium conflgurations 
of the system can be determined from these saddle-point equations. 

In current experiments, the temperature is well below the Fermi temperature Tp, and 
the system is often studied in the strong interaction regime where the magnitude of the 
scattering length can be much bigger than the interparticle distance. In this case, the 



interparticle distance becomes the most important length scale, and the thermodynamic 
properties of the Fermi gas are unitary [l9| . To the leading order, the free energy density of 
the normal state is given by 
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Ln = d^r -— CArns (r) + K(r)n(r) 
J 2m 



(5) 

where V is the trapping potential and c^ is a dimensionless constant. In the experiments ^, 
|20|, about 25 percent reduction of the energy by the interaction was observed, corresponding 
to Cat = 4.3. 

In the following, we study the phase-slip phenomenon in the unitary limit near the su- 
perfluid transition temperature, in the region which is likely to be accessible in experiments. 
The coefficients a and f3 in the unitary limit are functions of the particle density n and 
temperature T. Near the transition temperature, the coefficient a is approximately given 
by 

a = CaksiT -Tc)-, 

where Ca is a dimensionless constant, T is the temperature, and ks is the Boltzmann con- 
stant. In the unitary limit, the transition temperature is given by ksTc = ctJi n'^^^ /{2m), 
where ct is a dimensionless constant. The theoretical estimation of Tq is about 0.5 to 0.2 
iimes the Fermi temperature, and the corresponding constant ct is between 4.8 and 1.9 
2lL 122, |2^. The coefficients P in equation (J2I) is also a function of the fermion density, 
P = cph'^/[2m.n^^^], where c^g is a dimensionless constant. Currently there are no theoret- 



ical estimations of the constants Ca and C/j for the strong interaction in the unitary limit. 
However, the two constants are available for the case of weak interactions 2J], 

Ca = 1.474, (6) 

cp = 2.944. (7) 

The renormalization of these two constants by the strong interaction may not be significant, 
because the renormalization reduces the Fermi energy only by about 25 percent |^, |2i| . 

In the unitary limit near the transition temperature, the saddle point equations (0 EJ 
are given by 

- ^VV(r) + a(r)0(r) + /3(r)|</)(r)|V(r) = 0, (8) 
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where /x is the chemical potential. In equation (jH)), the term quartic in the order parameter 
has been dropped, because it is very small compared to other terms. In the ground state, 
the phase of the order parameter is same everywhere. Both the order parameter 0o(r) and 
the fermion density no{r) of the ground state can be determined from these equations. The 
order parameter of the ground state varies slowly. From equation (jH)), near the center of the 
trap, it is approximately given by 
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2mCakB{Tc - T)nj^{0) , 



where the origin is located at the center of the trap and 0o is chosen to be positive. 

In a phase-slip configuration, the phase of the order parameter changes by tt from one 
side of the trap to the other side. With this boundary condition, the order parameter 0i(r) 
and the fermion density ni(r) can also be obtained from the saddle-point equations (jHI) and 
0. Near the transition temperature, the density of the order parameter |0p is much less 
than the fermion density. From equation Q, the difference between the fermion density in 
the phase-slip configuration and that in ground state can be obtained, 

An(r) = ni(r) -no(r) ^ — — -An^{r), (11) 

OCtv 

where An^ is defined as An^{r) = |0i(r)p — |0o(r)p. 

In a quasi-lD trap, the system is tightly trapped in two directions, and the low-energy 
dynamics occurs in the other direction. Even in a trap which is not quasi-lD such as the 
„.e . Ref. Ud. as long as the widtK of the phase sl.p ,s „,ueh .ess than the t.ap s.e. the 
dynamics of the phase slip is approximately one dimensional near the center of the trap. In 
this case, equation (jH)) is reduced to a one- dimensional form, 

- ^^^ + ^i^M^) + /3(x)|0(x)|V(x) = 0, (12) 

and the order parameter of the phase-slip state has an analytic expression, 

0i(x) ?« 0o(O) tanh(Kx), (13) 



where x is the coordinate along the axial direction, k = J—ma{0)/h, and (f)i{x) is chosen to 
be real. As plotted in Fig. ((T)), in a phase-slip configuration, the order parameter vanishes 
in the middle. There is a phase difference of vr between the order parameter at one side and 



that at the other side. Far away from the middle, the magnitude of the order parameter 
approaches the value of the ground state. Note that the position of the phase slip does not 
change with time, because the state given by equation (fT^ does not carry any current. 

Near the center, according to equation (jllj) . the density change due to the phase slip is 
given by 

An(x) ^ -^^(f)l{0)sech\Kx), (14) 

OCn 

and in terms of the density and temperature, 

An{x) ^ ^^^^^^kB{T - Tc)nl^\0)sedi''{Kx). (15) 

The magnitude of the density change |An(a;)| reaches the maximum at the center and 
vanishes far away. The width of the phase slip in the real space is of the order of l/zt, where 



K = ^CamkB{Tc - T)/h. (16) 

In experiments, if k and the temperature can be measured, then the transition temperature 
Tc can be determined. Right at the center, the density change is proportional to the ratio 
of the condensation-energy density Cc to the normal-state-energy density Cn, 

An(0) 3TcCc{0) 



no{0) 5{T-Tc)Cn{0)^ 

which is equal to 



(17) 



An(0) _ 6cTclmkBiT-Tc) , . 

noiO) 5c,,c,h\ri0) ' ^ ^ 

where jCc{0) = a(O)|0o(O)P/2 and Cn{0) = cntiq (0). Thus the density profile of a phase 

slip can provide important information about the properties of the superfluid ground state 

such as the transition temperature and the condensation energy. 

In experiments, the method to detect the phase slip is the measurement of the density 

profile. However, the trap also causes the inhomogeneity of the density distribution. Both 

the phase slip and the trap lead to the quadratic dependence of the density on the distance 

away from the center, but they have different signs. Whether or not the phase slip shows up 

as a dip in the density profile is determined by the second derivative of the particle density 

at the center. The particle density can be written as the sum of the particle density of the 

ground-state and the density difference between the phase-slip state and the ground state, 

ni{x) = An{x) + nQ{x) . (19) 



The second derivative of the density change due to the phase shp is given by 

Near Tc, the density profile of the ground state is approximately given by the Thomas- Fermi 
approximation, 

—CNnUx) + V{x)-fx = 0, (21) 

and the second derivative of the density is approximately given by 

Qrn 1 Qm Ij) 1 

<(0) = -^nUO)r'{0) = -^n|(0), (22) 

on cn on cn 

where u^ is the axial frequency of the trap. Whether the phase slip appears as a dip in the 
density profile depends on the ratio of the two second derivatives, 



An"{0) 



n'o'(O) 



AcTclk%{T-Tc 



c 



(23) 



If this ratio is bigger than one, thenjhere is a dip in the density profile. 

In the experiment by Chin et al. [12], the excitation gap was observed at the Fermi energy 
Ep/ks = 3.2/iK with the axial frequency uJx/{2tt) = 24IIz. The transition temperature 
Tc was estimated to be about a quarter of the Fermi temperature |23], corresponding to 
ct = 2.4. Currently there are no better estimation of the constants c^ and C/3 than the values 
given by equations (jHI) and ((Zj). With these numbers, even if the temperature is lower than 
the transition temperature by just one percent, Tq — T = O.OlTc, the ratio of the two density 
derivatives given in equation ()23|) is already much larger than one, An"{0)/n'^{0) ^ 6.8 x 10^. 
Thus the phase slip can be identified by the dip at the center of the trap in the density profile 
for the vast range of the temperature below the transition temperature. In Fig. Q, the 
density profile of a phase slip at temperature T = 0.9Tc is plotted near the center of trap. 
At the center, the particle density of the phase slip is smaller than that of the ground state 
by about 12 percent, which is very likely to be detected in the experiment. Although better 
estimations of the constants such as c^ and c^ can improve these calculations, but it will 
probably not lead to any change greater than an order of a magnitude, as the interaction 
only changes the normal-state energy by about 25 percent in the unitary limit j^, |20|. The 
measurement of the constants Ca and Cf3 in the phase-slip experiment can also serve as a test 
for microscopic theories. 



The strong interaction is the reason for the unitary properties of a Fermi gas. However, 
when the Fermi gas is very close to a Feshbach resonance, there is a crossover from the 
superfluid state (BCS-type) to the molecular-condensate state (BEC). In the crossover re- 
gion, fermion atoms coexist with diatomic molecules, and properties of this mixture cannot 
be described by the unitary properties of the fermions alone, because of the interaction 
between atoms and molecules. Although in the crossover region many properties of the 
system are unclear, the phase slip can still be generated by the phase- imprinting method, 
because in both the superfluid and molecular-condensate states the phase symmetry of the 
wavef unction is broken. 

At the other side of the crossover, in the molecular-condensate state with a weak inter- 
action, the dynamics is well described by the Gross-Pitaevski equation (GP-equation) 

'^^ = (-|^V2 + ^)^ + 9m\i^\'^, (24) 

where ip is the wavefunction of the molecular condensate, gm is the coupling constant gm = 
h am/2m, and a^ is the scattering length of the molecules. The soliton solution of the GP- 
equation has been well studied. Near the center of the trap, the stationary soliton solution 
is approximately given by 

tpix, t) ^ ^0 tanh(K„x)e-2^A'tM^ (25) 

where the chemical potential is given by /i = gm\'ipo\'^/'2 and Hm = 2y/mfl/h. In this con- 
figuration, the molecular density vanishes at the center and the phase difference of the 
wave-function across the center is given by vr. The width of the phase slip is of the order 
of l/nm from which both the coupling constant gm and the scattering length am can be 
extracted. In Ref. J2a|, am was estimated to be about 0.6 times the scattering length of the 
atoms. The phase slip in the molecular-condensate state is much easier to be identified in 
the density profile because the density vanishes at the center. Comparing the properties of 
the phase slip at both sides of the BEC-BCS crossover, we expect that the particle density 
at the center drops dramatically as the system is tuned from the superfluid state to the 
molecular-condensate state, due to the increase of the molecular density. 

In conclusion, we predict that the phase slip can be generated in a superfluid Fermi gas 
by the phase imprinting technique. The presence of the phase slip can be detected by the 
measurement of the density profile. The phase slip can provide an excellent direct evidence 

8 



of the superfluid state, because phase coherence is one of the key signatures of superfluidity. 
The density profile of a phase slip also offers important information about the superfiuid 
state, such as the condensation energy and the transition temperature. We would like to 
acknowledge support from NSFC under Grant No. 90303008 and support from SRF for 
ROCS, SEM. 
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FIG. 1: The order parameter (j)i{x) of a phase-slip configuration near the center of trap. The 
width of the phase shp is much smaher than the trap size. The order parameter vanishes in the 
middle. Across the middle, the phase difference of the order parameter is given by vr. Away from 
the center, the magnitude of the order parameter approaches the value of the ground state. 
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FIG. 2: The density profile of a phase slip for the experimental conditions given in Ref. 1121 at 

the temperature T = 0.9Tc where Tq ~ 0.25rp, Tp = 3.4|uK, and ujx = 24Hz. The units of the 

1 
two axes are given by 1/k ^ l/[0.39no (0)] and no(0) ss 3.3 x 10^'^/cm'^. The value of no(0) is 

calculated from the Thomas-Fermi approximation. The width of the phase slip is of the order of 

1/k which is much smaller than the trap size. 
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